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We discuss the various ways in which primordial black holes may have formed in the 
early Universe and how the effects of such black holes can be used to place constraints on 
cosmological models. We show that such constraints may be severely modified if the value 
of the gravitational "constant" G varies with cosmological epoch, a possibility which arises 
in many scenarios for the early Universe. The nature of the modification depends upon 
whether the value of G near a black hole maintains the value it had at its formation epoch 
(corresponding to gravitational memory) or whether it tracks the background cosmological 
value. This is still uncertain but we discuss various approaches which might help to resolve 
the issue. 

> 

^NJ , §1. Introduction 

o 

£^ ■ It is well known that primordial black holes (PBHs) could have formed in the 

early Universe. 1 )' 2 ) A comparison of the cosmological density at any time after the 
Big Bang with the density associated with a black hole shows that PBHs would 
have of order the particle horizon mass at their formation epoch. PBHs could thus 
span an enormous mass range: those formed at the Planck time (10 43 s) would have 
the Planck mass (10 _5 g), whereas those formed at 1 s would be as large as 1O 5 M0, 
comparable to the mass of the holes thought to reside in galactic nuclei. 

PBHs would could arise in various ways. 3 ) Since the early Universe is unlikely 
to have been exactly Friedmann, they would form most naturally from initial inho- 
mogeneities but they might also form through other mechanisms at a cosmological 
phase transition. We discuss these various formation mechanisms in Section 2. PBHs 
could also have various types of cosmological consequence. Those larger than 10 15 g 
could contribute to the dark matter density and might even produce observable mi- 
crolensing effects. Those smaller than this would have evaporated by the Hawking 
mechanism 4 ^ and could contribute to the flux of cosmic rays. 

Studying the formation and cosmological consequences of PBHs is important 
because it enables one to place constraints on the early Universe (eg. on the spectrum 
of density fluctuations and the nature of any phase transitions) . Indeed PBHs serve 
as a probe of times much earlier than that associated with most other "relicts" of 
the Big Bang. While photons decoupled at 10 6 y, neutrinos at 1 s and WIMPs at 
lCU 10 s, non-evaporating PBHs go back to 10 _23 s and evaporating ones all the way 
back to the Planck time. Therefore even if PBHs never formed, their non-existence 
gives interesting information. We review and update these constraints in Section 2. 
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The main purpose of this paper is to examine how the PBH constraints are 
modified if the value of the gravitational "constant" G was different at early times. 
As reviewed in Section 3, this idea has a long history and should no longer be 
regarded as exotic. It arises in various scalar-tensor theories of gravity and these are 
a natural setting for many currently popular models of the early Universe. A number 
of astrophysical constraints suggest the G could not have varied much since the epoch 
of cosmological nucleosynthesis but PBHs would have formed much earlier than this. 
Homogeneous cosmological solutions in such theories have been much studied and 
generally lead to a simple power-law dependence of G on t. However, inhomogenous 
models (such as would be associated with PBH formation) would lead to a value of 
G varying both temporally and spatially and are much less well understood. 

Black hole formation and evaporation could be greatly modified in variable-G 
cosmologies, since many of their properties (eg. their radius and Hawking tempera- 
ture) depend explicitly on G. However, the nature of the modification depends upon 
the extent to which a PBH preserves the value of G at its formation epoch rather 
than following the background cosmological value. Barrow 5 ) first drew attention to 
this problem and introduced two possibilities: "scenario A", where G everywhere 
maintains the background value and "scenario B" , where the value of G stays con- 
stant near the black hole but evolves far away from it. This last possibility he termed 
"gravitational memory" . There would be interesting modifications to the cosmolog- 
ical consequences of PBH formation in both cases but they would be more dramatic 
in scenario B. 

Barrow & Carr 6 ) considered the implications of these two scenarios in detail and 
we will review some of their conclusions in Section 3. However, they did not address 
the issue of gravitational memory itself and we turn to this in Section 4. Although 
the problem is still not resolved, we will report on several interesting developments. 
These are described in more detail elsewhere 7 ** and exploit two important equiva- 
lences. The first equivalence is between scalar-tensor theories and general relativity 
with a scalar field. This comes about because these two theories are related by a 
conformal transformation and this means that the problem of gravitational memory 
can be probed by investigating the formation and evolution of a black hole in a 
cosmological background containing a scalar field. This leads to four variants of the 
gravitational memory scenario, which we describe in Section 5. The second equiva- 
lence is between a scalar field and a stiff fluid (with equation of state p = p). This 
allows one to relate the issue of gravitational memory to the problem of black hole 
accretion in a Universe containing a stiff fluid. This problem has already been inves- 
tigated by several people and, as discussed in Section 6, their findings lead to several 
new insights. 

§2. PBH formation and constraints on the early Universe 

One of the most important reasons for studying PBHs is that it enables one to 
place limits on the spectrum of density fluctuations in the early Universe. This is 
because, if the PBHs form directly from density perturbations, the fraction of regions 
undergoing collapse at any epoch is determined by the root-mean-square amplitude 
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e of the fluctuations entering the horizon at that epoch and the equation of state 
p = 7P (0 < 7 < 1). One usually expects a radiation equation of state (7 = 1/3) 
in the early Universe. In order to collapse against the pressure, an over dense region 
must be larger than the Jeans length at maximum expansion and this is just ^7 
times the horizon size. On the other hand, it cannot be larger than the horizon size, 
else it would form a separate closed universe and not be part of our Universe. 8 ) 

This has two important implications. Firstly, PBHs forming at time t should 
have of order the horizon mass then: 

M(t) « ^ « 10 15 ( I6 4 r -) g. (2-1) 

Secondly, for a region destined to collapse to a PBH, one requires the fractional 
overdensity at the horizon epoch, 5, to exceed 7. Providing the desnity fluctuations 
have a Gaussian distribution and are spherically symmetric, one can infer that the 
fraction of regions of mass M which collapse is 9 ) 



0{M) ~ e(M) exp 



7 2 



2e{Mf 



(2-2) 



where e(M) is the value of e when the horizon mass is M. The PBHs can have 
an extended mass spectrum only if the fluctuations are scale-invariant (i.e. with e 
independent of M) but this is expected in many scenarios. Recent hydrodynamical 
calculations for the 7 = 1/3 case have refined these estimates somewhat: Niemeyer 
& Jedamzik 10 ) find that one needs 5 > 0.8 rather than 5 > 0.3 to ensure PBH 
formation, and Shibata & Sasaki 11 ) reach similar conclusions, but the same basic 
picture still applies. 

Another interesting development has been the application of "critical phenom- 
ena" to PBH formation. Studies of the collapse of various types of spherically 
symmetric matter fields have shown that there is always a critical solution which 
separates those configurations which form a black hole from those which disperse 
to an asymptotically flat state. The configurations are described by some index p 
and, as the critical index p c is approached, the black hole mass is found to scale as 
(p — Pc) r ' for some exponent r/. This effect was first discovered for scalar fields 12 ) 
but subsequently demonstrated for radiation 13 ) and then more general fluids with 
equation of state p = jp. 14 )' 15 ) 

In all these studies the spacetime was assumed to be asymptotically flat. How- 
ever, Niemeyer &; Jedamzik 16 ) have recently applied the same idea to study black 
hole formation in asymptotically Friedmann models and have found similar results. 
For a variety of initial density perturbation profiles, they find that the relationship 
between the PBH mass and the the horizon-scale density perturbation has the form 

M = KMh(5 — <5 C ) 7 (2-3) 

where Mh is the horizon mass and the constants are in the range 0.34 < 7 < 0.37, 
2.4 < K < 11.9 and 0.67 < 5 C < 0.71 for the various configurations. Since M — ► as 
5 — > <5 C , this suggests that PBHs may be much smaller than the particle horizon at 
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formation (although it is clear that a fluid description must break down if they are 
too small) and it also modifies the mass spectrum. 17 )' 1S ) 

These developments might be regarded as refinements of the original scenario 
since they all presuppose a relationship of the form (2.2). However, in some situations 
eqn (2.2) would fail qualitatively. For example, PBHs would form more easily if the 
equation of state of the Universe were ever soft (7 <Cl). This might apply if there 
was a phase transition which channelled the mass of the Universe into non-relativistic 
particles or which temporally reduced the pressure. In this case, only those regions 
which are sufficiently spherically symmetric at maximum expansion can undergo 
collapse; the dependence of (3 on e would then be weaker than indicated by eqn (2.2) 
but there would still be a unique relationship between the two parameters. 19 ) 

The fluctuations required to make the PBHs may either be primordial or they 
may arise spontaneously at some epoch. One natural source of fluctuations would 
be inflation 20 ^' 21 ) and, in this context, e(M) depends implicitly on the inflationary 
potential. Indeed many people have studied PBH formation in inflationary scenarios 
as an important way of constraining this potential. 22 ) " 2S ) Recently the Gaussian 
assumption has been questioned in the inflationary context, 29 )' 30 ) so eqn (2.2) may 
not apply, but one still finds that (3 depends very sensitively on e. 

Some formation mechanisms for PBHs do not depend on having primordial fluc- 
tuations at all. For example, at any spontaneously broken symmetry epoch, PBHs 
might form through the collisions of bubbles of broken symmetry. 31 ) " 33 ) PBHs might 
also form spontaneously through the collapse of cosmic strings. 34 ) " 3S ) In these cases 
(3(M) depends, not on e(M), but on other cosmological parameters, such the bubble 
formation rate or the string mass-per-length. 

In all these scenarios, the current density parameter ftpBH associated with PBHs 
which form at a redshift z or time t is related to (3 by 9 ) 

-1/2 / M n-1/2 

^pbh = P(2r(1 + z)^ 10 6 /? (-J « 10 18 /? (j^r g ) (2-4) 

where Or ~ 10~ 4 is the density parameter of the microwave background and we have 
used eqn (2.1). The (1 + z) factor arises because the radiation density scales as (1 + 
z) 4 , whereas the PBH density scales as (1 + z) 3 . Any limit on £2pbh therefore places 
a constraint on (3(M) and the constraints are summarized in Fig. 1. The constraint 
for non-evaporating mass ranges above 10 15 g comes from requiring Qpbh < 1- Much 
stronger constraints are associated with PBHs which were smaller than this, since 
they would have evaporated by now. For example, the constraints below 10 6 g are 
based on the (not necessarily secure) assumption that evaporating PBHs leave stable 
Planck mass relics, in which case these relics are required to have less than the critical 
density. 23 )' 39 )' 40 ) Other constraints are associated with the generation of entropy, 
modifications to the cosmological production of light elements and the contribution 
to the cosmological 7-ray background. 

The constraints in Fig. 1 are discussed in detail by Carr et al. 23 ) but we note 
that Kohri & Yokoyama 41 ) have recently improved the constraints on /3(10 8 — 10 10 <7) 
which come from cosmological nucleosynthesis considerations. Here we wish to em- 
phasize that the strongest constraint is the cosmic ray limit associated with the 
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10 15 g PBHs evaporating at the present epoch. 42 ) " 46 ) Indeed the recent detection 
of a Galactic 7-ray background 47 ) , measurements of the antiproton flux, 48 ) and the 
discovery of very short period 7-ray burts, 49 ) may even provide positive evidence 
for such PBHs. This is discussed in detail elsewhere. 50 ) The constraints on (5{M) 
can be converted into constraints on e(M) using eqn (2.2) and these are shown in 
Fig. 2. Also shown here are the (non-PBH) constraints associated with the spec- 
tral distortions in the cosmic microwave background induced by the dissipation of 
intermediate scale density perturbations and the COBE quadrupole measurement, 
as well as lines corresponding to various slopes in the e(M) relationship. 

Finally it should be emphasized that there has been particular interest recently 
in whether PBHs could have formed at the quark-hadron phase transition at 10~ 5 s. 
This is because the horizon mass is of order IMq then, so such PBHs would naturally 
have the sort of mass required to explain the MACHO microlensing results. 51 - ) One 
might expect PBHs to form more easily at that epoch because of a temporary soft- 
ening of the equation of state. If the QCD phase transition is assumed to be of 1st 
order, then hydrodynamical calculations show that the value of 5 required for PBH 
formation is indeed reduced below the value which pertains in the radiation case. 52 ) 
This means that PBH formation will be strongly enhanced at the QCD epoch, with 
the mass distribution being peaked around the horizon mass then. 18 ) 

Since the focus of this meeting is gravitational waves, it should be emphasized 
that one of the interesting implications of the PBH MACHO scenario is the possible 
existence of a halo population of binary black holes. 53 ) In this case, there could be 
10 8 binaries inside 50 kpc and some of these could be coalescing due to gravitational 
radiation losses at the present epoch. 54 ) Current interferometers (LIGO, VIRGO, 
GEO, TAMA) could detect such coalescences within 50 Mpc, corresponding to a 
few events per year. Future space-borne interferometers (such as LISA or OMEGA) 
might detect 100 coalescences per year. If the associated gravitational waves were 
detected, it would provide a unique probe of the halo distribution (eg. its density 
profile and core radius) . 55 ) 

§3. Cosmology in varying-G theories 

Most variable-G scenarios associate the gravitational "constant" with some form 
of scalar field <p. This notion has its roots in Kaluza-Klein theory, in which a 
scalar field appears in the metric component 555 associated with the 5th dimension. 
Einstein-Maxwell theory then requires that this field be related to G. 56 ) Although 
this was assumed constant in the original Kaluza-Klein theory, Dirac 57 ) noted the 
the ratio of the electric to gravitational force between to protons (e 2 /Gto 2 ) and the 
ratio of the age of the Universe to the atomic timescale (t/t a ) and the square-root of 

the number of particles in the Universe (yj M/m p ) are all comparable and of order 

10 40 . This unlikely coincidence led him to propose that these relationships must 
always apply, which requires 



Goct -1 , GM/R~1, 
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Fig. 1. Constraints on /3(M) 



Fig. 2. Constraints on e(M) 

where R ~ ct is the horizon scale. The first condition led Jordan 58 ) to propose a 
theory in which the scalar field in Kaluza-Klein theory is a function of both space 
and time, and this then implies that G ~ has the same property. The second 
condition implies the Mach-type relationship <j) ~ M/R, which suggests 59 -* that <fi is a 
solution of the wave equation □</> ~ p. This motivated Brans-Dicke (BD) theory, 60 ^ 
in which the Einstein-Hilbert Lagrangian is replaced by 

L = 4>R — -J^^ u + L m , (3-2) 
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where L m is the matter Langrangian and the constant oj is the BD parameter. The 
potential 4> then satisfies 

where T is the trace of the matter stress-energy tensor, and this has the required 
Machian form. Since (f> must have a contribution from local sources of the form 
Ei(mi/ri), this entails a violation of the Strong Equivalence Principle. In order to 
test this, the PPN formalism was introduced. Applications of this test in a variety 
of astrophysical situations (involving the solar system, the binary pulsar and white 
dwarf cooling) currently require |cj| > 500, which implies that the deviations from 
general relativity can only ever be small in BD theory. 61 ) 

The introduction of generalized scalar-tensor theories, 62 ) ~ 64 ) in which oj is itself 
a function of (j), led to a considerably broader range of variable-G theories. In partic- 
ular, it permitted the possibility that oj may have been small at early times (allowing 
noticeable variations of G then) even if it is large today. In the last decade interest 
in such theories has been revitalized as a result of early Universe studies. Inflation 
theory 65 ) has made the introduction of scalar fields almost mandatory and extended 
inflation specifically requires a model in which G varies. 33 ) In higher dimensional 
Kaluza-Klein-type cosmologies, the variation in the sizes of the extra dimensions 
also naturally leads to a variation in G. 66 ^ " 68 ^ . The currently popular low energy 
string cosmologies necessarily involve a scalar (dilaton) field 69 ) and bosonic super- 
string theory, in particular, leads 70 ) to a Lagrangian of the form (3.2) with oj = — 1. 

The intimate connection between dilatons, inflatons and scalar-tensor theory 
arises because one can always transform from the (physical) Jordan frame to the 
Einstein frame, in which the Lagrangian has the standard Einstein-Hilbert form 71 ) 

L = R- 2V>Vv/r y + L m . (3-4) 

Here the new scalar field ip is defined by 

^={-^) J (3-5) 

and the barred (Einstein) metric and gravitational constant are related to the original 
(Jordan) ones by 

= A(cf>)%„, G=[l + a 2 (^]A(<p) 2 G, a = A'/A, (3-6) 

where the function A{4>) specifies a conformal transformation. Thus scalar-tensor 
theory can be interpreted as general relativity plus a scalar field. 

The behaviour of homogeneous cosmological models in BD theory is well under- 
stood. 72 ) Their crucial feature is that they are vacuum-dominated at early times but 
always tend towards the general relativistic solution during the radiation-dominated 
era. This is a consequence of the fact that the radiation energy-momentum tensor is 
trace-free [i.e. T = in eqn (3.3)]. This means that the full radiation solution can 
be approximated by joining a BD vacuum solution to a general relativistic radiation 
solution at some time t\, which may be regarded as a free parameter of the theory. 
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However, when the matter density becomes greater than the radiation density at 
t e ~ 10 11 s, the equation of state becomes that of dust (p = 0) and G begins to vary 
again. For a k = model, one can show that in the three eras 6 ) 

G = G (t /t e ) n , a«^ 3 (t>t e ) (3-7) 

G = G e = G {t /t e ) n , aoct 1 / 2 {h < t < t e ) (3-8) 

G = G e (Vil)~ {n+V4 " + "" )/2 , a OC t (2-n-V4n+n^)/6 (f < h) (3-9) 

where Go is the value of G at the current time to, n = 2/(4 + 3lo) and (to/t e ) ~ 10 6 . 

Since the BD coupling constant is constrained by \lo\ > 500, which implies 
|n| < 0.001, eqns (3.7) to (3.9) imply that the deviations from general relativity are 
never large if the value of n is always the same. However, for reasons explained 
below, it is also interesting to consider BD models in which n and lo can violate the 
current constraints. This allows considerably more exotic behaviour, especially in 
the vacuum-dominated era. In particular, we note the following features: 6 ) models 
with lo < —3/2 are probably excluded because the energy density of the scalar field 
is negative; models with —3/2 < lo < —4/3 undergo power-law inflation during the 
vacuum-dominated era because the exponent of t in the expression for a in eqn (3.9) 
exceeds 1; models with — 4/3 < lo < can bounce during the vacuum-dominated era 
because this exponent is negative. 

The behaviour of cosmological models in more general scalar-tensor theories 
depends on the form of lo(4>) but they still retain the feature that the general rela- 
tivistic solution is a late-time attractor during the radiation era. Since one requires 
G ~ Go to 10% at the epoch of primordial nucleosynthesis, 72 -* one needs the vacuum- 
dominated phase to end at some time t v < 1 s. The theory approaches general rela- 
tivity in the weak field limit only if lo — > oo and lo' /lo 3 — ► (where a prime denotes 
d/d(j)) but Lo(cp) is otherwise unconstrained. Barrow & Carr consider a toy model in 
which 

2lo + 3 = 2/3(1 - (t)/(t> c r a (3-10) 
where a and (3 are constants. This leads to 

2w + 3ocr a/(2 - a) , lo'/lo 3 oc t(i-2«)/(i-«) (t<t v ), (3-11) 

so one requires l/2<a<2in order to have lo — > oo and lo'/lo 3 — ► as t — > oo. In 
the q = 1 case, one finds 

G oct- 2X / (3 - x \ aott (1 - x V (3 - x \ A = ^3/(2/3) (*<*„). (3-12) 

During the vacuum-dominated era, such models can therefore be regarded as BD 
solutions in which lo is determined by the parameter (3 and unconstrained by any 
limits on lo at the present epoch. After t v , G is constant and one has the standard 
radiation-dominated or dust-dominated behaviour. 

On the assumption that a PBH of mass M has a temperature and mass-loss 
rate 4 ) 

T = (87TGM)- 1 , M w -(GM)- 2 , (3-13) 
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with G having the value G(t) in scenario A and G(M) in scenario B, Barrow & Carr 
calculate the evaporation time r for various values of the parameters n and t\ in BD 
theory. The results are shown in Fig. 3(a) for scenario A and Fig. 3(b) for scenario 
B. Here M* is the mass of a PBH evaporating at the present epoch, M e is the mass 
of a PBH evaporating at time t e and M cr u is the mass of a PBH evaporating at the 
present epoch in the standard (constant G) scenario. In scenario A with n < —1/2, 
there is a maximum mass of a PBH which can ever evaporate and this is denoted 
by Mqo. The results for the scalar-tensor with to((p) given by eqn (3.10) with a = 1 
are shown in Fig. 3(c) for scenario B with various values of the parameters A and 
t v . The corresponding modifications to the constraints on (3{M) in all three cases 
are shown in Fig. 3(d), which should be compared to Fig. 1. 



Fig. 3. Dependence of the PBH evaporation time r on initial mass M in (a) BD theory with scenario 
A, (b) BD theory with scenario B, (c) scalar-tensor theory with a = 1 and scenario B. Also 
shown are (d) the modifications to the constraints on f3(M) in these cases. 



§4. Black holes in scalar-tensor theory 

In BD theory Hawking 73 ) showed that, providing the weak energy condition 
holds, the gradient of <\> must be zero everywhere for stationary, asymptotically flat 
black holes. This means that such black holes are identical to those in general 



10 



B. Carr and C. Goymer 



relativity. This result can be generalized to all scalar-tensor theories and suggests 
that such theories are in agreement with the "no-hair" theorem. 

Numerical calculations support this theorem. 74 ) " 76 ) An Oppenhiemer-Snyder 
type collapse reveals outgoing scalar gravitational radiation, which radiates away 
the scalar mass. When all the scalar mass is lost, the black hole settles down with 
a constant scalar field to the Schwarzschild form. The scalar radiation generated 
in this way would be of great interest in itself. As shown by Harada et al., 76 ) if a 
black hole of mass M collapses, the ratio of the amplitudes of the scalar and tensor 
gravitational waves is 

~ io- 3 u— r 1 f~(—) 1/2 (^-) (4-i) 

h T \f)\ a J VT^' 1 \McJ \l5kmj' 14 ij 

where a is the initial size of the region, e is its eccentricity and / is the frequency. 
There might also be a stochastic background of scalar gravitational waves generated 
by the formation of PBHs in the early Universe, although its density would have 
been diluted by now because of redshift effects. 

It should be pointed out that the simulations of scalar collapse also show a rather 
unusual feature: 74 )' 75 ) the area of the event horizon decreases for a time, violating 
the area theorem, and it may also lie inside the apparent horizon during this period. 
This occurs for all values of u and is due to the fact that the scalar field violates 
the weak energy condition. This condition holds in the Einstein frame as long as 
oj > —3/2 but it need not be true in the physical Jordan frame. Examples have been 
found 77 ) " 79 ) of black holes in BD theory which have scalar hair but these arise only 
for certain ranges of u: u < —4/3 in 77 ), u < —3/2 in 78 ) and —5/2 < oj < —3/2 
in 79 ). It is unclear that these can represent the end state of realistic gravitational 
collapse. 

In addressing the question of which of Barrow's scenarios A or B is most plau- 
sible, it should be stressed that the scalar no hair theorem has only been proved for 
asymptotically flat spacetimes, so it is not clear that it also applies in the asymp- 
totically Friedmann case. While the no hair theorem suggests that <p should tend to 
a locally constant value (close to the black hole), it is not obvious that this needs to 
be the asymptotic cosmological value. Indeed, since the homogeneizing of <f> is only 
ensured by scalar wave emission, one might infer that this can only be achieved on 
scales less than the particle horizon. 

The only way to determine what happens is to seek a precise mathematical model 
for a black hole in a cosmological background. One approach is to try matching the 
black hole and cosmological solutions over some boundary S. An example of such a 
matching in general relativity is the Einstein-Straus or "Swiss cheese" model. 80 ) Here 
a Friedmann exterior is matched with a general spherically symmetric interior. If 
there is no scalar field, it turns out that the latter has to be the static Schwarzschild 
solution but the situation may be more complicated in the present context due to the 
presence of scalar gravitational radiation. In general one can show that the following 
continuity conditions must apply at U: 

M=0, [G^n< i n v ] = [G^uW] = 0, [<P] = 0, [cf>^] = (4-2) 
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where and u M are 4- vectors normal and tangent to £, respectively. 81 ) 

Unfortunately, it turns out that an Einstein-Straus type solution does not exist 
in BD theory. This is because the only way to satisfy the junction conditions (4.2) 
is if is spatially and temporally constant, which is just the general relativistic 
case. However, Oppenheimer-Snyder collapse has been investigated 76 )' 82 ) in which 
a ball of dust described by a k = +1 Friedmann interior and a Schwarzschild exterior 
collapses to a black hole. In these calculations the scalar field is taken to be constant 
before the collapse and its back-reaction on the metric is assumed to be always 
negligible. It is found that, as the collapse proceeds, a scalar gravitational wave 
propagates outwards before the scalar field settles down to being constant again. In 
principle, it should be possible to extend this approach to the present problem by 
attaching a k = +1 Friedman interior to a k = Friedman exterior. 

Jacobsen 83 ) has addressed the problem analytically by looking for a spherically 
symmetric solution which represents a black hole in an asymptotically Friedmann 
model in which satisfies the wave equation 

{-g tt d 2 t + -±=d r [^-g g rr d r ^j 0(t, r) = 0. (4-3) 

As r — > oo, he assumes that asymptotes to 

Mt) ~ 4>ct (4-4) 

in the Einstein frame. This presupposes that the black hole event horizon is much 
smaller than the particle horizon, so that the cosmological timescale is much longer 
than the black hole timescale, in which case C in eqn (4.4) can be regarded as 
constant. He then seeks a solution which is a combination of a homogeneous part 
01 (t) and a time-independent part 02 {r), where 

01 = 4> c t, 02 = ln(l - r H /r). (4-5) 

Both these terms diverge at the event horizon (r#) but one can find a combination 
of them which is regular there: 

03 = 01 + ^o020c = 4> c [v -r - r H ln(r/r H )} (4-6) 

where v is advanced "tortoise" time, together with a superposition of waves falling 
into the black hole or dispersing to infinity. The equipotential = 03 intersects the 
event horizon at v = vh and infinity at t = t^. Since v = v h and and t = too 
themselves intersect at r rj 1.5GM, he infers that there is little lag between the 
value of at the event horizon and particle horizon, which means that the memory 
can only be weak. 

Although Jacobsen's analysis demonstrates that a solution with little memory 
does exist, it must be emphasized that he has really put this feature into the solution 
at the outset by assuming the black hole is much smaller than the particle horizon. 
This assumption is inappropriate if the black hole has a size comparable to the particle 
horizon at formation (as applies for a PBH), so the field around a collapsing region 
need not necessarily evolve to the form (4.6). Another rather crucial feature of his 
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analysis is that he claims that, while the PBH mass is nearly constant in the Einstein 
frame, it increases as 

M = A(<t>y l M oc cft 1/2 (4-7) 

in the Jordan frame. If this were correct, it would invalidate the analysis of Barrow 
& Carr and hence Figs (3). However, we would argue that it is not clear in which 
frame the mass should be taken to be constant. 

§5. Variations of gravitational memory 

We now consider the question of what happens to black holes in BD or more 
general scalar-tensor theories during the evolution of the Universe. In general the 
present value of the scalar field, <p(to), will be different from its value when the black 
hole was formed, 4>(tf). In the following discussion, we will characterize the degree 
of gravitational memory by comparing the value of the scalar field at the black hole 
event horizon (4>eh) and the cosmological particle horizon {4>ph)- We first consider 
the two extreme situations described by Barrow 5 ) : 

Scenario A : 4>EH{t) = <ftpH(t) for all t (5 - l) 

A Schwarzschild black hole forms at time tf with its event horizon radius being 
Rf = 2G(tf)M. If G(t) evolves with time, then the black hole adjusts quasi-statically 
through a sequence of Schwarzschild states approximated by R = 2G(t)M, see Fig. 
4(a). In this scenario there are no stationary black holes when G(t) is changing and 
no gravitational memory. 

Scenario B : <j>EH(t) = <pEH{tf) for all t (5-2) 

A Schwarzschild black hole of size Rf forms at time tf and, while G(t) equals the 
evolving background value beyond some scale-length R m > Rf, it remains constant 
within R m , see Fig. 4(b). In this case the black hole size is determined by G{tf) 
even at the present epoch and this means that the region R < R m has a memory of 
the gravitational "constant" at the time of its formation. 

Neither of these scenarios can be completely realistic since they both assume 
that (ft is homogeneous almost everywhere. However, even if (ft were homogeneous 
initially, one would expect it to become inhomogeneous as collapse proceeds. Indeed, 
if the background value is increasing (as usually applies), one would expect (ft in 
the collapsing region to become first larger than the background value on a local 
dynamical timescale and then smaller than it on a cosmological timescale. Such 
behaviour would necessarily entail a variation of (ft in space as well as time. Solutions 
have been found for specific to in which (ft is inhomogeneous 84 \ so such models are 
clearly possible in principle. We must also allow for the possibility that (ft may vary 
interior to R m but on a slower or faster timescale than the background. We therefore 
propose two further scenarios: 
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Scenario C : 



\4>EH{t)\>\4>PH{t)\ for all t 



(5-3) 



where the dot represents a time derivative. This implies that the scalar field evolves 
faster at the event horizon than at the particle horizon until it eventually becomes 
homogeneous, see Fig. 4(c). We describe this as short-term gravitational memory 
and it reduces to scenario A as the timescale to become homogeneous tends to zero. 
This would apply, for example, if 4> were to change on the dynamical timescales of 
the black hole since this is usually less than that of the cosmological background. 



This implies that (j) evolves faster at the particle horizon than the event horizon, 
see Fig. 4(d). We describe this is as weak gravitational memory and it reduces to 
scenario B when the left-hand-side of eqn (5-4) is zero. In this case, the evolution of 
(j) is again dominated by the black hole inside some length-scale R m . Note that, in 
either this scenario or the last one, the length-scale R m need not be fixed, since it 
could either grow or shrink as scalar gravitational radiation propagates. A particu- 
lar example of this, to which we return shortly, would be self-similar gravitational 
memory, in which the ratio of <f>EH to 4>pu always remains the same. 



Scenario D : 



\4>EH(t)\ <\(/> PH (t)\ for all t 



(5-4) 



Fig. 4. Different possiible forms for the evolution of the scalar field profile <j>(r) for (a) no memory, 
(b) strong memory, (c) short-term memory, (d) weak memory. 
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If <f>, or equivalently G, is inhomogeneous, then one has to consider carefully 
what is meant by the Hawking temperature. In general relativity the temperature of 
a Schwarzschild black hole is given by eqn (3.13). However, in a scalar-tensor theory 
G = G(t,r), so one must decide which value of G is appropriate. Since the particle 
creation responsible for the effect takes place near the event horizon, perhaps the 
local value Geh should be used, giving 

T H = (SitGehM)- 1 . (5-5) 

On the other hand, since the radiation is only measured asymptotically, one might 
argue that the asymptotic value Goo is appropriate, giving 

T H = (SttGooM)- 1 . (5-6) 

It is not certain whether eqn (5-5) or (5-6) should be used. If eqn (5-5) applies, then 
the temperature depends crucially on whether or not gravitational memory exists. 
However, if eqn (5-6) applies, the temperature is unaffected by this and just depends 
on the background. 

Gravitational memory has also been investigated in the context of boson stars 85 ^ . 
A boson star is the analogue of a neutron star, formed when a large collection of 
bosonic particles become gravitationally bound 86 )> 87 ). In this case, weak gravita- 
tional memory is already known to occur. An interesting feature of this is that two 
boson stars of the same mass may differ in other physical properties (eg. radius), 
depending on when they formed. However, the situation may be very different for a 
black hole since this has an event horizon. 



§6. Gravitational memory and the accretion of a stiff fluid 

In general relativity there is an equivalence between a scalar field and a stiff 
fluid and we now show how this can be exploited in studying gravitational memory. 
In the Einstein frame, the energy momentum tensor for a perfect fluid is 

T^u = (p + p)u^u u + g^p (6-1) 

where is the velocity of the fluid. If we define a velocity field by 

4>ix 



this gives 



(6-2) 



Tr = J£±g&. +I *.. (6-3) 
By comparing this to the energy-momentum tensor for a scalar field, we find that 

p = p = -\g pa ^ P ^, (6-4) 

so we have a stiff fluid. This equivalence applies provided that the derivative of 
the scalar field is timelike. Otherwise the velocity field defined in (6-2) would be 
imaginary. 
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This is relevant to the gravitational memory problem because we can now inter- 
pret the various scenarios discussed in Section 5 in terms of the accretion of a stiff 
fluid. If the black hole does not accrete at all or accretes very little, this will cor- 
respond to strong or weak gravitational memory (scenarios B and D, respectively). 
However, if enough accretion occurs to homogenize eft, this will correspond to short- 
term gravitational memory (scenario C). The faster the accretion, the shorter the 
memory, so scenario A corresponds to the idealization in which homogenization is 
instantaneous. 

A simple Newtonian treatment ^ for a general fluid suggests that the accretion 
rate in the Einstein frame should be 

M = 4irpR 2 A v s , (6-5) 

where Ra = GM/v 2 is the accretion radius and v s is the sound-speed in the accreted 
fluid. For a stiff fluid, v s = c and R A = GM/c 2 , while p ~ l/(Gt 2 ) in a Friedmann 
universe at early times, so we have 

dM GM 2 
~dF ~ ~cW 

This can be integrated to give 

M w 

where Mf is the black hole mass at the time tf when it formed. If we assume that 
Mf = rjtf with rj < 1, then eqn (6-7) implies 

M -> Mfil-r])- 1 as t-»oo. (6-8) 

If rj <C 1, the black hole could not grow very much. However, if 77 is close to 1, which 
must be the case if v s ~ c, then the black hole could grow significantly. In particular, 
in the limit 77 = 1, eqn (6-7) implies M ~ t, so the black hole grows at the same rate 
as the universe. This simple calculation suggests that a black hole surrounded by a 
stiff fluid can accrete enough to grow at the same rate as the Universe. 

Since the above calculation neglects the effects of the cosmological expansion, 
one needs a relativistic calculation to check this. The Newtonian result suggests 
that one should look for a spherically symmetric self-similar solution, in which every 
dimensionless variable is a function of z = r/t, so that it is unchanged by the 
transformation t — ► at, r — > ar for any constant a. This problem has an interesting 
but rather convolved history. By looking for a black hole solution attached to an 
exact Friedmann solution via a sonic point, Carr & Hawking first showed that there 
is no such solution for a radiation fluid 8 ) and the argument can be extended to a 
general p = jp fluid with < 7 < 1. Lin et al. 87 ) subsequently claimed that there 
is such a solution in the special case 7 = 1. However, Bicknell & Henriksen 88 ) then 
showed that this solution is unphysical, in that the density gradient diverges at the 
event horizon. The solution can be completed only by attaching a Vaidya ingoing 



(6-6) 



(6-7) 
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radiation solution interior to some surface (i.e. the scalar field has to turn into a 
null fluid). If one regards this as unphysical, this suggests that the black hole must 
soon become much smaller than the particle horizon, after which eqn (6.7) implies 
there will be very little further accretion. Therefore the stiff fluid analysis suggests 
that there should be at least weak gravitational memory. 
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